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Abstract 

Based on the trace anomaly for the energy-momentum tensor, an effective 
theory for the thermodynamics of the deconfining phase, and by assuming 
the asymptotic behavior to be determined by one-loop perturbation theory 
we compute the nonperturbative beta function for the fundamental coupling 
g in SU(2) and SU(3) Yang-Mills theory. With increasing temperature we 
observe a very rapid approach to the perturbative running. The Landau pole 
is nonperturbatively screened. 



1 Introduction 



Knowledge on the resolution dependence (running) of the fundamental coupling in 
Yang-Mills theories was first obtained within a perturbative setting pQ. This is im- 
portant because the perturbative renormalizability of these theories [2] , see also [3] , 
states that (apart from a wave-function renormalization) the coupling remains the 
only parameter of the theory in the process of successive perturbative elimination of 
quantum fluctuations. The proof of the perturbative renormalizability of Yang-Mills 
theories belongs to the deepest and most far-reaching theoretical results of the last 
century. As a consequence, the notion of asymptotic freedom was established for the 
strong interactions governed by the gauge group SU(3) pQ. The asymptotic freedom 
of Quantum Chromodynamics allows for a controlled small- coupling expansion of 
correlation functions around the conformal limit which takes place at an infinitely 
large resolution. 

Nonperturbative investigations on the running of the gauge coupling have been 
carried out in the framework of the exact renormalization group equation (ERGE) 
and the approach via the Dyson-Schwinger equations, see [U E] and references 
therein. As it seems, these results indicate that the perturbatively derived Landau 
pole is regularized by nonperturbative effects. The purpose of the present article is 
to re-investigate this issue in an independent way. 

On the one hand, we use a nonperturbative approach to SU(2) and SU(3) Yang- 
Mills thermodynamics [B] , which predicts for the deconfining phase the existence of 
a unique, maximal resolution in terms of the modulus |0|(T) of an emergent, adjoint 
Higgs fielcQ- On the other hand, we combine it with a nonperturbative definition 
of the gauge-coupling evolution via the trace anomaly for the energy-momentum 
tensor at finite temperature. This allows for an extraction of the nonperturbative 
rate of change of the fundamental coupling g under a variation of temperature: The 
beta-function. 

As a boundary condition we require that the high-temperature behavior of the 
beta function is in accord with the perturbative situation. The so-extracted law 
governing the running of the coupling is in agreement with that obtained in one-loop 
perturbation theory (when setting the renormalization scale equal to temperature) 
except closely above the critical temperature. 

The article is organized as follows. In Sec. [2] we briefly review the trace anomaly 
for the energy-momentum tensor and discuss its validity at finite temperature. Based 
on the trace anomaly and an effective theory for the thermodynamics of the decon- 
fining phase the derivation of the nonperturbative evolution equation for the fun- 
damental coupling in dependence of temperature is performed for SU(2) and SU(3) 
pure Yang-Mills theories in Sec. [33 We discuss the high-temperature limit to make 
contact with perturbation theory and subsequently solve the evolution equations. 
In Sec. H] we summarize and discuss our results. 

1 The spatially infinitely extended thermalized Yang-Mills system can be considered as an ex- 
tended vertex which selfconsistently picks its own scale of maximal resolution [BJ. 
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2 Trace anomaly 



The trace anomaly for the energy- momentum tensor 9^, which is considered an 
operator identity, reads [7J [8] 



= Pjil p a F a (1) 



where = d^A a v — d v A a — g f abc A b ^A c v is the field strength appearing in the fun- 
damental Lagrangian of the Yang-Mills theory, Cyai = ~\(F^ U ) 2 , the f abc are the 
structure constants of the Lie algebra, and (3 is given by the right-hand side of the 
evolution equation for the fundamental gauge coupling g: 

^d,g = (3(g). (2) 

In Eq. (J2J) the mass scale \x refers to the resolution that is applied to the process at 
which the strength of the coupling g is extracted. In contrast to the chiral anomaly 
[7j, which is not renormalized because of its topological nature, the trace anomaly 
exhibits two resolution- dependent factors: The /3-function divided by g and the 
average of F^F^. 

When solving [3(g) = fid^g = —bg 3 (b = valid at one-loop order and zero- 
temperature) a Landau pole \x = Al, which roughly is identified with the Yang-Mills 
scale A (A L ~ A), occurs: 

A L = no exp (--Lj . (3) 
Eq. ([3]) then implies the well-known form for the running coupling: 

9 = l + 261n(^) = 2bHfT) ■ (4) 

Let us now turn to the case of finite temperature: We work in a flat Euclidean 
spacetime with time r constrained as < r < 1/T. In j9j the one-loop zero- 
temperature expression in Eq. (j3J) was used to argue for the validitjH of Eq. (JTJ when 
taking a thermal average. 

Performing a thermal average over Eq. (Q]), we obtain 

P - = ^ (F; u F a ^) T , (5) 



2 The grand canonical potential density, y = — P, must be expressible as A 4 f(A/T) where / is 
a dimensionless function of its dimensionless argument. Taking the derivative of y with respect to 
the exponent in Eq. ^ , the thermal average over both sides of Eq. fl]) appears with the one-loop 
level expression for (3. The authors of [3] then claim that there is no difficulty to extend this result 
to an arbitrary order in perturbation theory. Here we will consider the one-loop situation only. 
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where p and p describe the energy density and the pressure of the thermalized Yang- 
Mills system. Notice that in Eq. (JS]) two scales enter: the temperature T, at which 
the thermal averages are calculated, and the scale /i, which is the resolution at which 
(3(g) = pd^g is evaluated [T0l[TT]. In our approach [6], [12], [13] the scales p and T are 
not independent but functionally related. 



3 Nonperturbative /3-function 

3.1 Effective SU(2) Yang-Mills thermodynamics 

We now turn to the effective theory for SU(2)-thermodynamics in the deconfining 
phase. There are topologically trivial, coarse-grained gauge fields a M entering in 
the effective field strength G a ^ v = d^a") — d u (a a ^) — ee abc a b ^al (to be distinguished 
from the fundamental field strength F^ v in Cym = — i(-^) 2 )' an d there is an inert, 
adjoint scalar field 0, which together with a pure-gauge configuration represents 
the thermal ground state emerging from a spatial average over interacting calorons 
and anticalorons of topological charge modulus \Q\ = 1. The effective coupling e 
is temperature dependent (to be distinguished form the fundamental coupling g). 
The dependence e = e(T) follows from thermodynamical selfconsistency, see below. 
The effective Lagrangian for the description of SU (2)- Yang-Mills thermodynamics 
in the deconfining phase (T > T c = A c A/27r, A c = 13.87) and in unitary gauge reads 

mm- 

£tt s = \ (GTW) 2 + MT? \<P\ 2 (K 1 )) 2 + «)) 2 ) + . (6) 

The modulus of the adjoint scalar field \4>\ depends on the Yang-Mills scale A and 

on temperature T as \<p\ = \J ^f,- The length is the minimal length down 

to which the thermalized system looks spatially homogeneous. In other words, the 
spatial average over interacting calorons and anticalorons selfconsistently saturates 
at this length scale. The quantum fluctuations ' are massive in a temperature 
dependent way, m 2 = m(T) 2 = 4e 2 |0| 2 , while the gauge mode ajP , stays massless 
(dynamical gauge symmetry breaking: SU(2) — > U(l)). 
We work with the following dimensionless quantities 

- P - P , 2ttT ... m(T) e(T).,. 

where p and p are the energy density and the pressure due to the Lagrangian (jBJ), 
and the function a = a(X) is introduced for later use. 

The energy density and pressure p and p are a sum of three contributions 

P = P3+ Pl,2 + Pgs, P = P3+ Pl,2 + Pgs , (8) 



3 



(1 2) 

where the subscript 1,2 is understood as a sum over the two massive modes , the 

(3) 

subscript 3 refers to the massless mode ajj, , and the subscript gs labels the ground- 
state contribution. When expressing p and p as functions of the dimensionless 
temperature A, one obtains at one-loop (accurate on the 0.1%-level [6| [14"]): 



7T 
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3 f°° . x 2 Vx 2 + a 2 _ 2(2tt) 



2 q /"OC 

2 90' ^ = ~V 2 / x2rfxln l 1 - e~ vOT ) , P gs = - j5 flS • (10) 
Imposing the validity of the thermodynamical Legendre transformation 



r^dP dp , . 

and substituting the expressions f[9"j)- f[TUj) into ffTTl) . we arrive at the following differ- 
ential equation for a = a(A): 

6A 3 / , da \ . , , . 

A— + a ) aD(a), (12) 



(2tt) 6 V rfA 
f°° x 2 1 

B(0) = I dX ^ + a ^^'-l ' a(Xm) = - <13) 

For a sufficiently large initial value Aj n the solution for a(A) is independent on Aj n : 
a low-temperature attractor with a logarithmic pole at A c = 13.87 is seen to exist. 
The effective coupling is given as e = e(A) = a(A)A 3 / 2 /47r, and exhibits a plateau 
e = v8~7r for A A c . In fact, 

■w = ^ (") 

is a solution of the differential equation ( Tl2|) for a <C 1, that is, for A ^> A c . For 
plots and the discussion of the thermodynamical quantities we refer to [6], and for 
the discussion of the linear growth of {0^) T to [T5] . 

3.2 Nonperturbative running coupling g(T): SU (2)-c&se 

We now use the effective Lagrangian ([6]) to evaluate the two relevant elements of 
the trace-anomaly equation 



P-3P=^-(F°„F^) T . (15) 



Namely, 

(i) p — 3p is evaluated from Eqs. (l9l- (TT0l . 

(ii) The expectation value (F^F a,Atiy ) T is the average action density in euclidean 
space (a possible definition of the gluon condensate as discussed in [T3~| [16]). 
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We evaluate the average action density by utilizing the effective Lagrangian (J6]) 
considering that the part of the fundamental field strength F® u , which enters the 
ground-state physics described by the effective theorj|§, suffers a wave-function renor- 
malization. We have 

(£ YM ) T = ~ (F^F a ^) T = f{g) (Clt B ) T = f{g)p gs (16) 

with p gs = X 4 p gs = in A 3 T . The function f(g) will be fixed by requiring that for 
A 3> A c the fundamental coupling g runs in agreement with perturbation theory. 
Within the effective theory we can relate the natural (not externally imposed) 

resolution scale p to temperature T. We have p = \<f>\ = \J [HI [121 H3]- That 
is, the thermalized Yang-Mills system acts like a spatially extended vertex being 
probed with a selfconsistently adjusting resolution \<p\ as soon as a temperature T 
is provided. Fluctuations, that would be resolved at p > \<f>\, are integrated out in 
the effective theory. Thus we have: 

13(g) =pd,g = -2Td T g = -2{3 T (g) , (17) 

where /3r{g) = Tdrg- Notice that (3{g) = pd^g is a positive quantity: In fact, 
the trace anomaly p — 3p and (F® u F a, ' JjU ') are both positive. The fact that the 
resolution p = \<fi\ decreases for increasing T generates a negative Pr(g) m accord 
with asymptotic freedom (T ^> T c ). 

Taking into Eqs. ffT5l) . ffT6l) , and ffTTj) . we have 

MA) - ^ = -^f( 9 ). (18) 
^pgs g 

The function h(X), see also [T5] . is plotted in Fig.[TJ Notice that h approaches the 
value | for A ^> A c thus implying that (0^) T = p — 3p = 6p gs = 247rA 3 T for T ^> T C) . 
This simple high-T behavior allows to determine the function f(g) analytically. We 
require that the perturbative result Pr(g) = — bg 3 {b = N = 2) holds for g <C 1 
(or T ^> T c ). Then Eq. (TT5|) implies that 




(19) 

The evolution equation (fIS|) for g = g(\) is recast as: 

M9) = ~\h HX)g 3 ^d x g = - 2 -b ^-g 3 . (20) 



3 On the one hand, it is straight-forward to show that on-shell excitations do not contribute to 
the thermal average (£dcc- c ff) T in Eq. © on the one-loop level. On the other hand, in the effective 
theory quantum fluctuations make a negligible contribution as compared to that of the ground 
state [BJ. 
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Figure 1: The function h(X), defined in Eq. ( 1T81) . plotted for the SU(2) (gray curve) 
and for the SU(3) (black curve) Yang-Mills theories. 

From the behavior of h(X) we can immediately infer two interesting properties: 

a) The function h(X) ~ | for A > 5A C . That is, the perturbative equation Prig) = 
— bg 3 is valid all the way down to 5A C . The range of validity of the perturbative 
treatment for the determination of g(X) is thus even larger than one naively would 
think. Indeed, as we will see later, the perturbative solution is very similar to the 
nonperturbative one even down to temperatures ~ 1.2 T c . 

b) The function h(X) slowly decreases for decreasing temperatures thus effectively 
lowering the coefficient b in the perturbative beta function. Therefore a mild screen- 
ing of the perturbative Landau pole occurs. 

To solve Eq. fTSOj) an initial condition must be specified. As discussed in the 
previous subsection, the effective coupling constant e(A) shows a logarithmic pole 
at the critical temperature A c = 13.87. We therefore impose that the fundamental 
coupling g(X) also diverges at the deconfining temperature. In this way the solution 
g = g(X) is uniquely fixed. Fig. [2] shows this solution and the perturbative solution 
gp(X) of Eq. PJ (by setting /i = A). The boundary condition for gp(X) is determined 
by imposing that g and gp coincide at large values of A. In practice, we can set g(X = 
10 A c ) = gp(X = 10 A c ) (another choice of the matching at high T clearly would lead 
to very similar results). Notice that in Fig. [2] the perturbative Landau pole is to the 
right of A c . Also, the nonperturbative coupling g is screened as compared to gp: 
g(X) is always below gp(X). Before moving on to the SU(3) case one comment is 
in order. Namely, our determination of the wave-function renormalization f(g) is 
based on the one-loop leading, asymptotic behavior of the beta function. Higher- 
order perturbative corrections, however, depend on the adopted renormalization 
scheme. Demanding scheme-invariance of the nonperturbative coupling, we are left 
with the one-loop expression for f(g). 
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Figure 2: The running of the nonperturbative g(X) (undotted) and that of the 
perturbartive counterpart gp(X) (dotted) for the SU(2)case. 



3.3 Nonperturbative running coupling g(T): SU (3)-case 

Here the effective thermodynamic description follows the same lines as in the SU(2) 
case, see [6]. We only report on some relevant formulas and briefly discuss their 
consequences. The modulus of the scalar field <j) is exactly the same. As shown 
in [6] out of the eight coarse-grained gauge modes four acquire a mass mi = e\<p\ 
(contributing to p and p by p\ and pi), two a mass m 2 = 2e \4>\ {p2 and P2), and two 
stay massless (^3 and ps). Explicitly, we have (a = mi/T): 
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Thermo dynamical selfconsistency, Eq. (flTT) . implies that 

12A 3 



1 



(2vr) f 



+ a] (aD(a) + 2aD(2a)) . 
dX J 



(23) 
(24) 

(25) 



8 ^2 



7r"A 3 / 2 , and the effective cou- 



The asymptotic solution to Eq. (1251) reads a(A) 
pling reaches a plateau value of e = for A ^> A c . The value for A c , where the 
effective coupling e possesses a logarithmic pole, is A c = 9.475 [6]. The asymptotic 
value of the function h(X), defined in (ITS]) , is 3/2 just like in the SU(2) case. The 
function h(X) is plotted in Fig. [[J and an analogous behavior to the SU(2) curve is 
evident. 
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Figure 3: The running of the nonperturbative coupling g(X) (undotted) and the 
perturbartive counterpart gp(X) (dotted) in the SU(3) case. 

Equations (fT9l) and (120]) hold for b = N = 3. The qualitative discussion is 
very similar to the SU(2) case. The nonperturbative and the perturbative solutions 
are plotted in Fig.0 where a screening of the perturbative Landau pole is observed. 
The similarity of the SU(2) and SU(3) cases is remarkable. 

4 Summary and Conclusions 

We have computed the beta-function for the fundamental coupling in SU(2) and 
SU(3) Yang-Mills theory by appealing to the trace anomaly for the energy-momentum 
tensor and an effective theory for the thermodynamics of the deconfining phase. The 
latter involves a adjoint scalar field <fi which emerges upon a coarse-graining process 
performed over interacting calorons and anticalorons [6] and together with a pure- 
gauge configuration represent the thermal ground state for the system. In contrast 
to earlier (perturbative) investigations, where a separate dependence of the coupling 
on the resolution \i and temperature T was employed [TUl [EE] , the effective theory 
dictates the relevant ii at a given T in terms of the modulus of the scalar field \<p\. 
Therefore T and ii are no longer independent variables, but they are functionally 



related: fi = \<f>\ = y where A is the Yang- Mills scale. The beta-function, 
defined by the rate of change of the fundamental coupling g when varying the reso- 
lution scale (a, thus also relates to the rate of change (3t of g when varying T. This 
enables a comparison with the one-loop, zero-temperature perturbative prediction 
(obtained by setting ii — T). In doing so we have assumed that for asymptoti- 
cally large temperatures the full beta function exhibits the universal, perturbative 
one-loop behavior. 

Remarkably, already for temperatures slightly greater than the critical temper- 
ature, we observe the behavior of (3t as predicted by one-loop perturbation theory 
[1]. Deviations become apparent only for smaller temperature, say below 1.5 T c . The 
nonperturbative effects generate a screening of the perturbative Landau pole: The 
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nonperturbative pole is slightly shifted to the left of the perturbative Landau pole 
for both SU(2) and SU(3). 

The fact that in our study the perturbative solution is valid up to temperatures 
slightly above T c relies on the properties of the function h(X) = as plotted 

in Fig.[T] (A = y" an d p gs is the explicit contribution of the ground state to the 
energy density). At high T the asymptotic behavior is h(X) ~ 3/2, and the per- 
turbative evolution equation is recovered. In the context of the present paper, the 
very presence of the nonperturbative ground state is ultimately responsible of this 
simple asymptotic behavior of h(X). Notice that the asymptotic behavior for h(X) 
implies that p — 3p = 6p gs = 24ttA 3 T. Thus the trace anomaly grows linearly with 
temperature for sufficiently high T, see also [TBJ [T71 [T8"] . 
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